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Abstract 

The non-perturbative nilpotent exact BRST symmetry of the Gribov-Zwanziger action in the Lan¬ 
dau gauge constructed in [T] is generalized to the case of Euclidean Yang-Mills theories quantized 
in the maximal Abelian gauge. The resulting diagonal gluon propagator is evaluating in dimensions 
D = 4,3, 2. In D = 4,3 a decoupling type behaviour is found in the infrared region, while in D = 2 a 
scaling type behaviour emerges. 


1 Introduction 


Since the seminal work [2], the issue of the Gribov copies has become a powerful tool in order to inves¬ 
tigate the behaviour of confining Yang-Mills theories in the non-perturbative infrared region. The deep 
progress [21 El SI El El El [H] done in the last four decades on the Gribov problem0in the Landau gauge, has 
resulted in the so-called Gribov-Zwanziger framework mmm and its refined version mmm, from 
which a local and renormalizable non-perturbative action taking into account the existence of Gribov 
copies has been constructed. So far, this set up has provided a large number of applications which have 
covered several aspects of confining Yang-Mills theories, namely: study of the gluon and ghost correlation 
functions mummmmm, investigation of the spectrum of the glueballs [33 IS], thermodynamics 
and phase transitions nzmamaEDiEiiEzzicai, supersymmetric theories Eliza, study of the confine¬ 
ment/deconfinement transition when Higgs fields are present [26J. 

Parallel to these developments, an intense and rich discussion on the important aspect of the relationship 
between the BRST symmetry and the Gribov problem has taken place, see [2H1 HJ ESI EH EZ EZ El 
IMl [Ml E71 [Ml [Ml IMl [Ml [Ml [M IM IM] for an overview. Needless to say, the BRST symmetry is a 
fundamental ingredient of the Faddeev-Popov quantization, allowing for a perturbative all order proof 
of the renormalizability of Yang-Mills theories as well as for the identification of the asymptotic Fock 
sub-spaces on which the S'-matrisd is unitarity. Unravelling the role of the BRST symmetry in the case 
of confining theories is believed to be a pivotal issue in order to understand the mechanism responsible 
for color confinement. 

*caprimarcio@gmail.com 
t diegofiorentinia@gmail.com 

* silvio. sorella@gmail. com 

x For a pedagogical introduction to the Gribov problem, see [9l 1101 . 

2 We refer here to non-confining Yang-Mills theories for which the elementary fields admit asymptotic states, so that a 
perturbative construction of the Ymatrix can be worked out. 
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Recently, the authors [Tj have been able to construct a non-perturbative nilpotent extension of the 
standard BRST operator which turns out to be an exact symmetry of the Gribov-Zwanziger action in 
Landau gauge. This is an important step towards a comprehension of the role of this symmetry at the 
non-perturbative level. The construction of the non-perturbative BRST operator outlined in pQ| has 
allowed for a nice and geometrical resolution of the Gribov problem in the class of the linear covariant 
gauges, with the important outcome that the correlation functions of the gauge invariant colorless opera¬ 
tors are independent of the gauge parameter, a feature also shared by the so-called Gribov parameter 7 2 . 
These are non-trivial results, given the well known difficulties of addressing the Gribov issue in the class 
of the linear covariant gauges [[46, 47], due to the lack of hermiticity of the Faddeev-Popov operator. 

In this work, we extend the construction of [lj to the case of the maximal Abelian gauge, which plays 
a central role in the dual superconductivity mechanism for confinement [481 4j9[ [50] • According to the 
dual superconductivity picture, QCD at low energies should behave as an Abelian theory in presence of 
monopoles, a feature referred as to Abelian dominance [5H [52j [53[ [Ml ESI [56l EH [58l ESI ESI EH E2 E3 
f6il 65]. The condensation of the monopoles would give rise to the formation of flux tubes which confine 
quarks (48j[49,, 50] . 

Concerning the Gribov problem in the maximal Abelian gauge, although the situation cannot be compared 
to that of the Landau gauge [HU El EiE], a few results are already available, see [Ml EH ESI ES [701 EH [72] , 
where the analogous of Zwanziger’s horizon function as well as of the Gribov-Zwanziger action and of its 
refined version have been constructed. A study of the maximal Abelian gauge within the context of the 
Schwinger-Dyson equations can be found in m- 

The paper is organized as follows. In Sect.2 we give a short summary of the construction of the non- 
perturbative BRST symmetry in the case of the Gribov-Zwanziger action in the Landau gauge. Sect.3 
is devoted to the generalization to the maximal Abelian gauge. In Sect.4 we discuss the refined action 
and its non-perturbative BRST symmetry. In Sect.5 we analyse the resulting diagonal gluon propagator 
in D = 4,3, 2 dimensions pointing out that, similarly to what happens in the case of the Landau gauge, 
while a decoupling behaviour takes place in D = 4,3, a scaling type behaviour emerges in D = 2. In 
Sect.6 we provide a detailed account of how the non-perturbative BRST symmetry and associated ac¬ 
tion can be cast in local form by introducing a suitable set of localizing fields. In Sect.7 we collect our 
conclusion. 


2 Brief summary of the Gribov-Zwanziger framework and of its non- 
perturbative BRST symmetry in Landau gauge 


For the benefit of the reader, let us give here a short summary of the Gribov-Zwanziger set up for 
SU(N ) Euclidean Yang-Mills theory in the Landau gauge, = 0, and of its non-perturbative BRST 

symmetry. This framework El 12 El El El [TO] enables us to implement the restriction in the functional 
integral to the Gribov region 0, defined as 

n = ( K\ d >‘ A l = o. Mab ( A ) > o} > (i) 

where M. ab is the Faddeev-Popov operator of the Landau gauge, i.e. 

M ab (A) = -5 ab d 2 + gf abc A c ^d^ with <9 M A“ = 0. (2) 

The restriction to the region 0 takes into account the existence of the Gribov copies which affect the 
Faddeev-Popov quantization scheme El 12 El El El ED]. According to E1I2E1E], for the partition function 
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of quantized Yang-Mills action in Landau gauge, one writes 


Z = [ [DA] 5(dA) det(Ad) e" SYM . 

JVL 


(3) 


The restriction of the domain of integration to the region hi can be effectively implemented by adding to 
the starting action an additional non-local term H(A), known as the horizon function. More precisely 

mmmm 


where 



[DA] 6(dA ) det(M) e" 5 ™ 


= J [DA] S(dA) det(.M) e -(s™+hH(A)-4vA(N*-i)) 
H(A) = g 2 J d A xd A y f abc A b ^x) [M” 1 (x, y)} ad f dec A^y), 


(4) 

(5) 


with [Ad -1 ] denoting the inverse of the Faddeev-Popov operator, eq.©. The mass parameter q 2 ap¬ 
pearing in expression Q is known as the Gribov parameter. It is determined in a self-consistent way by 
the gap equation mmmm 


(H) = 4V(N 2 - 1), (6) 

where the vacuum expectation values ( H ) has to be evaluated with the measure defined in eq. (J4]) ; V 
denotes the space-time volume. Expression ([J]) can be cast in a more suitable form by introducing a set 
of commuting ((j), <f>) and anticommuting (uj,ui) auxiliary fields [4j [5j E|, namely 

[ [DA] 5{dA ) det(M) e“ SYM = [ [D<4>] wh(N*-i))^ ( 7 ) 

.hi J 

where <h refers to all fields present and Sqz stands for the Gribov-Zwanziger action, i.e. 

Sgz = S FP + j d 4 x fa c [M(A)] ab ^ - Co; c [M(A)] ab u b ; + 91 2 f abc A*^ + $=)) , ( 8 ) 

with S F p denoting the Faddeev-Popov action of the Landau gauge 

Sfp = S YM + J d 4 x (b a d^A“ + <?dvD ab c b ) . ( 9 ) 

Notice that, in the expression ([9]) above, the field b a denotes the Lagrange multiplier enforcing the Landau 
gauge, while c a and c a are the Faddeev-Popov ghost and anti-ghost fields, respectively. 

In the local formulation, the gap equation (0 can be rewritten as 

|% = 0 , e~ V£v = J [L><f>] e-iSaz-^HN*-!)) , (10) 

where £ v denotes the vacuum energy. 
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In order to construct a non-perturbative BRST symmetry for the action (| 8 j). we follow [1J and intro¬ 
duce the non-local transverse gauge invariant fielqj A b d y A b = 0 , [75) 76, j77] 


Ah _ p 

-rt-fl — r LLV 


A v - ig 
d„ 


DA 
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+ 0(A 3 ) 
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, dA 


+ 0(A 3 ), 


( 11 ) 


where = [S^y — is the transverse projector. Expression dill) is left invariant by infinitesimal 

gauge transformations order by order, namely 


8Ap — 0 , 6Ap — —dX + ig[Afj,, A] . ( 12 ) 

Looking now at eq. (HU). one sees that a divergence dA is present in all higher order terms. Therefore, 
we can rewrite Zwanziger’s horizon function H(A) in terms of the invariant field A h as 

H(A) = H(A h ) - R{A){dA) (13) 


where R{A){dA) is a short-hand notation, R(A)(dA) = f d 4 xd 4 yR a (x,g)(dA a ) y , R(A) being an infinite 
non-local power series of A^. Thus, for the Gribov-Zwanziger action, we may write 


Sgz = S YM + J d 4 x [b a d u Al + Ad.Dfc^ + 7 4 H(A) 

= S YM + jd 4 x(b a d^Al + Ad^c b )+ 1 A H(A h )- 1 4 R{A){dA) 

= Sym + j d 4 x (b ah d„A“ + <fd^D ab c c ) + 7 4 P(AL h ) , (14) 

where the new field b ah stands for 

b h = b- r y 4 R(A) . (15) 

The use of the field b h enables us to write down an exact nilpotent non-perturbative BRST symmetry. 
Rewriting the Gribov-Zwanziger action by using the auxiliary fields ((/>, 4>,uj,Co), i.e. 

Sgz = Sym + J d 4 x (b ah d y A“ + Ad y D ab c b ) 

+ J d 4 x ^{M{A h )} ab ^ - u“[M{A h )] ab tf + gi 2 f abc Af{^ + <^ c )) , (16) 


it turns out that expression d! 6 D is left invariant by the nilpotent non-perturbative BRST transformation 

m 


S 7 2C“ 


~D ab c b , 

uah 


s^c a = g -r b A b c c , 

,ah 


s 7 2 b an = 0 


,ab 


^7 , = 0 > 

< b + 7 2 gf kpb A h ' k \M~\A h ) 


- pa 


3 We employ here a matrix notation, A^. = A^,' a T a , where {T a |. a = 1, ...(N 2 
SU(N), [T a ,T b ] = if abc T c . 


= 0 > (17) 

1), denote the hermitian generators of 
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with 


« 7 2=0, s^Sgz = 0 . (18) 

The operator s 7 2 depends explicitly on the non-perturbative Gribov parameter y 2 . As such, it is a non- 
perturbative extension of the usual BRST operator of the Faddeev-Popov theory, to which it reduces in 
the limit y 2 —>■ 0. 

Before discussing the generalization of this construction to the case of the maximal Abelian gauge, it is 
worth spending a few words on the nature of the non-local variable of expression (ED- As already 
pointed out, the field is left invariant order by order by infinitesimal gauge transformations. As a 
consequence, it is invariant under the action of the BRST operator s 7 2 , i.e. 

= 0 • ( 19 ) 

Equation (1191) might give the impression that the gauge field could create an invariant physical state 
carrying a color index. Nevertheless, in the present case, this possibility is ruled out by the presence of 
the Gribov parameter y 2 , which acts as a confining parameter for gluons in the non-perturbative infrared 
region. In fact, a look at the tree-level two-point correlation function (A^(k)A^(—k))GZ stemming from 
the Gribov-Zwanziger action ([16D learns that 

( A ,i{k)A„(-k)} G z = k4 + ^2 iV7 4 (V - ) ■ ( 20 ) 

Expression (1201) displays complex poles, a feature which prevents a particle interpretation. In other words, 
even if the field A^ is left invariant by the operator s 7 2 , it cannot be associated to a physical state due to 
the confining character of the Gribov horizon, encoded in the Gribov mass parameter y 2 . These consider¬ 
ations can be easily extended to the case of the refined-Gribov-Zwanziger theory Dunams] as well as to 
the maximal Abelian gauge. Of course, the field A 1 * could be introduced also in the case of a non-confining 
theory, like the electroweak theory SU(2) x 17(1), see [26]. In this case, the use of the variable A 1 * leads 
to a gauge invariant description of the physical vector bosons (IT + , W~, Z°), as discussed in details in m- 

The introduction of the field A^ and of the operator s 7 2 allows us to clarify better the physical meaning 
of the Gribov parameter y 2 . As it is apparent from expressions (jl4D . (1161) . the horizon function y 4 H(A h ) 
acquires now the meaning of a quantity invariant under the BRST operator s 7 2 . This is an interesting 
feature, enabling us to state that the Gribov parameter y 2 is not akin to an unphysical gauge parameter. 
In fact, from 

s 7 2 H{A h ) = 0 , (21) 

it immediately follows that 

( 22 ) 

meaning that y 2 is a physical parameter of the theory. As such, it will enter in the expressions of the 
gauge invariant correlation functions, deeply modifying their behavior in the infrared region, as explicitly 
reported, for example, in the evaluation of the spectral densities for the glueball spectrum nans]. 


3 Generalization to the maximal Abelian gauge 


In order to generalize the previous set up to the maximal Abelian gauge, let us start by fixing the 
notation and by reminding a few properties of this gauge. 
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3.1 The maximal Abelian gauge fixing and the horizon function 


We shall consider, for simplicity, the gauge group 517(2), whose generators T a , (o = 1,2,3), 


rj-iCL rj-tC 


= ie abc T c 


( 23 ) 


are chosen to be hermitian and to obey the orthonormality condition Tr (T a T b ) = 5 ab . Following 
SHI ESI [79|, we decompose the gauge field A M into off-diagonal and diagonal components 

A ,„ = ApT a = A«T a + A^T :3 , (24) 

where a = 1,2 and T 3 is the diagonal generator of the Cartan subgroup of 517(2). Analogously, decom¬ 
posing the held strength, we obtain 


rp _ T?a rpa _ jp'G. rpcx _i_ jp nn 3 

r fiv — r /JLV 1 — r /us 1 ' r in' 1 •) 

with the off-diagonal and diagonal components given, respectively, by 

F%, = DfAP-DfA?, 


/IV 


= d^A v - d v A^ ge^A^Al, 


where the covariant derivative Dpf is dehned with respect to the diagonal component A fi 

Df = - ge^Ap , e ap = . 

Also, the Yang-Mills action 

Sym = J d A x ^ (F^F^ + F ixv F ilv ) , 

is left invariant by the gauge transformations 

= -Dfuf-ge^Afr, 

SAa = -dnU - g£ af3 A a u)P . 


(25) 

(26) 

(27) 

(28) 


(29) 


The maximal Abelian gauge condition amounts to impose that the off-diagonal components A" of the 
gauge held obey the following non-linear condition 


D°-PAi = 0 . 


(30) 


Moreover, as it is apparent from the presence of the covariant derivative Dff, equation m allows 
for a residual local 17(1) invariance corresponding to the diagonal subgroup of 51/(2). This additional 
invariance has to be hxed by means of a further gauge condition on the diagonal component A /t , which 
is usually chosen to be of the Landau type, namely 


d^ = 0 


(31) 


The Faddeev-Popov operator, A4"Y corresponding to the gauge condition (1301) is easily derived, being 
given by 

= -DpDf - . (32) 

In order to construct the Faddeev-Popov action corresponding to the gauge conditions (Eqd,® , we 
proceed by introducing the standard nilpotent BRST transformations 


= -(DfP + ge^A^c), sA^ = -(5 M c + ge ap A«P) 


sc a = ge a P(Pc, 


sc = -e a Pc a P , 

2 




= b° 


sc = 


sb a = sb = 0 , 


(33) 


6 




where (c“,c, c“,c) are the Faddeev-Popov ghosts and ( b a ,b ) are the Lagrange multipliers implementing 
the gauge conditions (IHUD . (1HTT) . Further, we introduce the s-exact gauge fixing term 


Smag = s j d A x (?*D“ 0 A 0 + cd^A^j 

= f d A x (b a DfA 0 - c?*M a0 c 0 + ge^iD&Afo + bd^A,,. + cd^ (d^c + ge a0 A°c 0 ^ , 

(34) 

where A4 ab stands for the Faddeev-Popov operator (1321) . It is apparent that the Faddeev-Popov action 
of the maximal Abelian gauge 

Smag = Sym + Smag , (35) 

turns out to be BRST invariant, i.e. 

s Smag = 0 • (36) 

As any other covariant gauge, also the maximal Abelian gauge is plagued by the existence of Gribov 
copies, see refs. [80], j81j [82] for explicit examples of zero modes of the Faddeev-Popov operator (l32lh By 
restricting the integration in the functional integral to the region where the Faddeev-Popov operator AA a0 
is strictly positive, i.e. M. a0 > 0, a large number of copies could be eliminated, as proven in ESU23. 
Furthermore, in complete analogy with the case of the Landau gauge, this restriction is implemented by 
adding to the original Faddeev-Popov action, ea. ([35l) . a non-local horizon term which, in the case of the 
maximal Abelian gauge, turns out to be given by the expression [Ml EH E2HS3 

H M ag(A) = g 2 J d A x d 4 y A^{x)e a0 (Ad" 1 )" 5 (x,y)e 50 A^y). (37) 

Therefore, for the analogous of the Gribov-Zwanziger action in the maximal Abelian gauge, we have 

SaFag = Smag + 7 4 Hmag(A) , (38) 


where y 2 stands for the Gribov parameter of the maximal Abelian gauge. Proceeding as in the case of 
the Landau gauge, expression (|38[) can be cast in local form by introducing a pair of auxiliary bosonic 
fields, (4> < j' l 0 and a pair of auxiliary fermonic fields, {ugf, w^), namely 

S G m\g = SIFag + J d A x (tfM****? - + 91 2 e a0 (f# - 4>f) A,) . (39) 

As shown in [661 EZl [69UT0], the action S^ag ena bles us to implement the restriction in the functional 
integral to the Gribov region Qmag °f the maximal Abelian gauge, defined as 


^MAG = | , Aft, 


<9 uA^ = 0, DfA 0 = 0, M a0 (A) 


(DfD 50 +g >e a °E 0& A°X) >0 } - (40) 


Although the understanding of the Gribov issue in the maximal Abelian gauge cannot yet be compared 
to that reached in the Landau gauge, a few properties of the region £Imag have been already obtained. 
In particular, in [TO], it has been established that &mag is unbounded along the diagonal directions in 
field space. This feature seems to be consistent with the aforementioned Abelian dominance hypothesis, 
according to which the diagonal configurations, corresponding to the Abelian Cartan subgroup, should be 
the dominant configurations in the infrared. Moreover, in im it has been shown that when an Abelian 
configuration is gauge-transformed to the Landau gauge, it is mapped into a point of the boundary of the 
Gribov region 17 of the Landau gauge, ea. pUl) . i.e. into a point of the Gribov horizorfl These features 
give further support to the restriction of the domain of integration to the region Qmag- 


We can now address the issue of the existence of a nilpotent non-perturbative BRST symmetry for 
the action (1391) . 

4 See Sect.V of [73] . 
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3.2 The non-perturbative BRST symmetry for the maximal Abelian gauge 


In order to construct a non-perturbative BRST symmetry for the action (1391) . we proceed as in the 
case of the Landau gauge. Making use of the invariant field A j ). eq. (HU), we can rewrite expression (EU 
as 


Hmag(A) = H MAG (A h )-T(A)dA-T a (A)dA a 

= H MAG (A h ) - (-d^(A) + ge a ^ a (A)A^) A, - F a {A)DfA^ (41) 

where, as in the previous section, J r {A)dA, F a (A)dA a stand for a short-hand notation, i.e. J-(A)(dA ) = 
/ d 4 xd 4 yJ 7 (x,y)(dA) y and F a {A){dA a ) = f d 4 xd 4 yJ ra (x,y)(dA a ) y . Introducing as before the redefined 
Lagrange multipliers b h ,b h ' a 

b h = b- 1 4 T{A) + 1 4 j X dy^[g£^F a {A)A^ y 

b h ’ a = b a — q 4 J 7 " (A) , (42) 

we can rewrite the action (1381) as 

Sm Z ag = &YM + SMAc(b h , b h,a ) + 7 4 H M AG(A h ) , (43) 

with 

S MAG (b h ,b h ’ a ) = f d 4 x (b h ’ a DfAP - c a M a Pc p + ge^iD^A^c + b h d^ + cd„ (d^c + ge^A^Jj 

(44) 

Notice that equations (|42l) correspond to a change of variables in the functional integral in the 6-sector 
of the theory, the corresponding Jacobian being the unity. 

Introducing thus the auxiliary fileds (4>,4>,Co,oj), we obtain 

S G m Z ag = Sym + J d 4 x (t?M aS (A h )<f> s J > - QfM aS {A h )Jf + gi 2 e(ff - tff) A\ )’ 3 ) 

+ f d 4 x (b h ’ a DfAP - (?M a0 c p + ge^ip^A^c + b h d^ + cd, M (d^c + ge af} . 

(45) 

As in the case of the Landau gauge, expression (Hhl) enjoys a non-perturbative exact nilpotent BRST 
symmetry, namely 



= + ge^A^c) , s^A^ = 

-{dpC + ge^A^) 

S 7 2C“ 

= ge^c^c, s^c = ~£ a Pc a c? , 


S 7 2C° 

= b h,a , s^2C = b h , 


Sl 2b h ’ a 

= s 7 2 6 = 0 , 



= , s 7 = 0 , 



= 4>f + gi 2 ^A h ^(M-\A h )) 5a , 

s 7 2$f = 0 , 


and 

Sj2Sy2 = 0 , S 7 2 Sm Z ag = 0 • (47) 

Let us end this section by mentioning that, although we have considered the gauge group 51/(2), the 
above construction can be straightforwardly generalized to the gauge group SU ( N ) by making use of the 
results of m- 


4 The Refined Gribov-Zwanziger action in the maximal Abelian gauge 


As already mentioned, the Gribov-Zwanziger action in the Landau gauge, eq. csd, generalizes to 
its refined version nn nans]. The origin of the refined version stems for the observation that non- 
perturbative condensates of dimensions two, (A^' a A^,’ a ) and (<^V2 6 — w“ 6 w“ b ), are dynamically generated 
for non-vanishing Gribov parameter y 2 . The effective action taking into account the existence of these 
condensates is called the refined Gribov-Zwanziger action mmm- So far, the prediction of the refined 
action on the infrared behaviour of the gluon and ghost propagators are in very good agreement with the 
most recent lattice data, see [83, .54]. 


A refined version of the action (14511 in the maximal Abelian gauge has also been constructed in [69], 
where the analogous of the dimension two condensates of the Landau gauge have been introduced. More 
precisely, in the case of the maximal Abelian gauge, due to the splitting of the gauge field into diagonal 
and off-diagonal components, we have the following dimension two condensates: (A^ ,3 A^ 3 ) 

and u/jf), playing a different role at the dynamical level. 


In particular, the condensate (A^A^’ 01 ) provides a dynamical Yukawa type mass for the off-diagonal 
components of the gluon field jSSfiSB]- This gives support to the Abelian dominance hypothesis, accord¬ 
ing to which the off-diagonal components should decouple at very low energy scales [581 ISlffl 


On the other hand, the condensates (A^' 3 A^' 3 ) and {(fff <j)ff — ujfujf) determine the infrared behaviour 
of the diagonal gluon propagator. Following (69) . for the refined action which takes into account the 
dimension two condensates in the maximal Abelian gauge, we have 


qRGZ _ q gz , 
MAG ~ a MAG + 


m 


off 


/ 


d 4 xA h ’ a A h ’ a + 

r A*- 


d A x 





) + 


m 


diag A h ’ 3 A h ’ 3 


(49) 


with S AI Z AG given in expression (1451) . 


The parameters (m 2 ®, /r 2 , m(j iag ) have a dynamical origin, encoding the presence of the dimension two 

condensates (A^’“A^’ a ), (A 1 ^’ 3 A 1 ^’ 3 ) and ((jfjf (ffjf — ujjf wjf). The non-perturbative BRST nilpotent sym¬ 
metry introduced in the previous section generalizes almost immediately to the case of the refined action 
m ■ The only difference with respect to equations (146K is given by the transformation of the field , 
which now reads 


= 


+ 91 


2 £ 60^h,3 


M(A h ) + fj . 2 1 


l -1 


dot 


(50) 


Again, one checks that 

s 7 2s 7 2 = 0 , s 7 2*Smag = 0 • 


(51) 


In view of the analysis of the diagonal gluon propagator, it is instructive to present the explicit first-order 
evaluation of the condensates (A^’ 3 A^’ 3 ) and <j>^ — ujjf ufjf). 


5 As a consequence of the existence of the condensate (A^’ a A ^’ a ), the tree level off-diagonal propagator takes the form 

<A“(fc)A£(-fc)> = ^ ’ ( 48 ) 

where the dynamical mass parameter m% s is related to (A^' a A^' a ). Numerical lattice simulations [Sit 1521 [531154115511561 
f57l l58l l59l l60l El; E21 l63l 164 , 65] have given evidence that the off-diagonal mass m][ ff is large enough so as the off-diagonal 
gluon components are more suppressed than the corresponding diagonal components at very low energy scales. 
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At the first-order only the quadratic part of the action ([4511 is needed, namely 


(^mag) quad 


= J d A x Qa"(-<5^<9 2 + + ^A^-S^d 2 + d^d„)A^j 

+ J d 4 x (6“<9 M A“ + c“<9 2 c" + bd^A^ + cd 2 c ) 

+ J d 4 x - of(-d 2 )uf + gi 2 e^ (if - if) A,) . 


( 52 ) 


Further, we introduce the operators f d 4 xAf Af and J d 4 x(if if — ujfujf) in the action by coupling 
them to two constant sources J and a, and we define the vacuum functional £(a, J ) defined by 


a - V£{<t,J) _ 


[D§] c ~( S mag)<i™a- J / b$*)-o f d 4 x 


(53) 


where we used the fact that, at the first-order, 


J d 4 x AfAf = J d 4 x A (V - A v . (54) 

It is easy to check that the condensates {Af Af) and (if if — ufuf) are obtained by differentiating 
£(<t, J) with respect to the sources (J, a), which are set to zero at the end, i.e. 



d£(a, J) 
dJ 

d£(a, J ) 
da 


J=a =0 


J=a =0 


Employing dimensional regularization, a direct computation shows that 


£{a, J) 


(D- 1 ) 
2 


d D k 

(2ir) D 


In 



4 yV 
k 2 + J 



+ 


terms independ. from (J, a). 


(55) 


(56) 


Ea. ([55|) and ([56]) give thus 


(4>f4>f ~ ufujf) Blst _ oldei = - 2 yV(£> - 1 ) 


d D k 1 1 

(27t) d k 2 ( k 4 + 4g 2 y 4 ) 


(57) 


and 


i Ah,3 Ah,3\ 

X-*- 1 //, / first-order 


-4ff 2 7 4 (^ - 1) 


1 1 

(2ir) D k 2 (fc 4 + 4g 2 y 4 ) 


(58) 


Eg. (1571) and ea. (|58|) show that, already at first-order, both condensates (AfAf) and (if if — wfuf) 
are non-vanishing and proportional to the Gribov parameter y 2 . Notice also that both integrals in 
eas. (l571) . (]58[) are perfectly convergent in the ultraviolet region by power counting. We see thus that, in 
perfect analogy with the case of the Landau gauge mmm, dimension two condensates are automat¬ 
ically generated by the presence of the Gribov parameter y 2 . 
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Having introduced the refined action, eo. (149l) . it is helpful to add a few remarks on its meaning. To 
that aim, we integrate out the auxiliary fields , <tyf Luff), obtaining 


qRGZ 
°MAG 


+ 


Sym + <?V [ d A x d A y A h ^{x)e afi f ^M{A h ) + yu 2 l] ^ e s/ 3 A h /(y) 

J \ L / x,y 

f d*xA h ^A h ^ + f d 4 x -^^) + ^A h ^Af 


+ J d A x (b h ’ a DfAd - g*M afi cP + ge afi c a (D^A^)c + b h d^ + cd M (d^c + ge af} , 


(59) 


from which one sees that the starting horizon function (1371) gets replaced by the expression 


Hmag(A h 


2 4 

9 7 


/ 


A h 


M(A h ) + y 2 1 


i -i 


A h 


( 60 ) 


Nevertheless, one also observes that the Faddeev-Popov operator Ai a d(A) entering the Faddeev-Popov 
ghost sector of expression (|59D is left unchanged, being given by the term c a M a d(A)c@. Thus, the 
question which naturally arises is: how the properties of the Faddeev-Popov ghost correlation function 
(c a cd) = (Ad _1 ) a ^ are changed when the gluon sector has been modified according to ea. (l60D ? To answer 
this question we shall partly refer to the case of the Landau gauge, where the same situation is found. 
First, we underline that the parameter y, 2 appearing in the refined action (1491) is not free. In fact, from 
expression m , one realizes that the dimension two condensate (4>(f> — lulu) is proportional to the Gribov 
parameter y 2 . This means that the mass parameter y 2 has a dynamical origin, being generated by the 
Gribov parameter itself. Furthermore, as shown in great detail in |12] in the case of the Landau gauge, 
where an explicit check of the no-pole Gribov condition has been worked out, with the inclusion of the 
dimension two condensate (<fr(j) — LULu) the Faddeev-Popov ghost correlation function (cc) = (A4 -1 ) retains 
the fundamental property of being positive. This means that with the inclusion of the condensate one 
remains inside the Gribov region Q, i.e. the Gribov horizon is never crossed. Thus, in presence of the 
condensate, the Faddeev-Popov ghost stays positive, as required by the restriction to the Gribov region. 
Though, its infrared behaviour is now deeply changed. More precisely, with the inclusion of dimension 
two condensate it is no more enhanced in the infrared, behaving as (cc) 1/A: 2 , as opposed to the 

/c ~0 

enhanced behaviour 1/A: 4 observed in the absence of condensate. Needless to say, a ghost behaviour of 
the kind 1/A; 2 is in very good agreement with the most recent numerical lattice simulations, as well as 
with the studies based on the Schwinger-Dyson approach. Finally, we underline that the presence of the 
dimension two condensates turns out to be energetically favoured, as shown in m through the explicit 
evaluation of the effective potential. 

We are now ready to evaluate the tree level diagonal gluon propagator. This will be the topic of the next 
section. 


5 The diagonal gluon propagator in D = 4, 3, 2 dimensions 


The four-dimensional diagonal gluon propagator can be read off from the refined action (1491) . being 
given by 


{A^{k)A u {-k)) D=A 



k^ky \ _ k 2 + y 2 _ 

k 2 J k 4 + (m^ iag + y 2 )k 2 + y 2 m^ iag + 4 ff 2 y 4 ’ 


( 61 ) 
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where the parameters (/r 2 , m;(j iag ) 


are related to the dimension two condensates 


/* 2 ~ (ifr/ - G?u?) , >»L ~ (A™ A™) . 


(62) 


A few remarks are in order here. We notice that expression (16111 shares a deep similarity with the gluon 
propagator of the Landau gauge mmm- As in the Landau case, it turns out to be suppressed in the 
low momentum region, attaining a non-vanishing value at zero momentum, i.e. at k 2 = 0. This kind of 
behaviour is referred to as the decoupling solution. Moreover, expression (1611) is in good agreement with 
the lattice numerical simulations of the diagonal gluon propagator in momentum space reported in [641165] . 


Although not explicitly mentioned, the actions (1451) and (1491) retain their validity in arbitrary dimen¬ 
sion D = 4, 3, 2, so that we can access directly the diagonal gluon propagator in both D = 3,2. 


Looking at the expressions of the condensates in eas. (l56l) . (1571) . we see that they exist in D = 3, the 
corresponding integrals being perfectly convergent in both IR and UV regions. As a consequence, for the 
diagonal gluon propagator in D = 3 we get a decoupling type behaviour as well, i.e. 


{A^k)A v (-k)) D= 3 



k,,k v \ k 2 + fj 2 

k 2 ) k A + (m 2 iag + g 2 )k 2 + g 2 m 2 iag + 4 g 2 A ' 


(63) 


Furthermore, from expressions (1571) . (1551) . one realizes that the dimension two condensates (A 1 ^’ 3 A 1 ^’ 3 ), 
<j)^f — uifuf) cannot be safely introduced in D = 2 dimensions, due to the presence of infrared 
singularities, as it is apparent from the presence of the term l/k 2 in the integrand of expressions ()57[) . (1581) . 
Instead, in D = 2 a scaling behaviour given by a Gribov type propagator is exhibited by the diagonal 
gluons, namely 

(A„(k)A„(-k)) D . 2 = (V - LL) k4 , (64) 

which, unlike expressions (1611) . (1631) . vanishes at zero momentum. Again, this feature is in full agreement 
with what observed in Landau gauge, where a scaling type solution emerges in D = 2 dimensions 
IS7HBSUH]. The emergence of a scaling type solution in D = 2 for the diagonal propagator might have 
also strong consequences on the validity of the Abelian dominance hypothesis in D = 2, as recently 
advocated in [ 90] , 


6 Localization of the non-local field and of the BRST operator s 2 

In the previous sections use has been made of non-local expressions, i.e. both gauge field A jj, eq. m , 
and BRST operator s 7 2 , ea. (]46l) . are written in terms of non-local quantities. It is thus relevant to show 
here how they can be cast in a local form, so that the standard tools of the algebraic renormalization 
m can be employed to analyse the structure of the theory and of its symmetry content. 

Let us first focus on the non-local field A^ and on the non-local part of the action (1451) . which we 
rewrite as 

Sm Z ag = s ym + S~mag 

+ f Ax ( b h ' Q DfA <1 - + ge a P c a (D^ 5 Afyc + b h d^ + c<9 M (d^c + g£ af 3 A°c^ . 

(65) 
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where Sj^ AG stands for the non-local expression 


07 _ 

°MAG ~ 


d 4 x (4>fM a5 {A h )^ - 0 jfM a \A h )ujff + 91 2 e^ A' 


h, 3 

p 


( 66 ) 


The field A^ can be localized by introducing an auxiliary Stueckelberg field £“ [92], namely 

= c%r a = h!(AiT a )h + -h%h, 

h = e^“ T \ 


(67) 


where £ a is a local dimensionless Stueckelberg field. The local field C^ can be now expanded in terms of 
the Stueckelberg field £“, yielding 


Cl = Al- - 9 -f abc (?D c X + 0(f). 


( 68 ) 


> m s 2 J s ^ 

Following [92], for the BRST transformation of h, one gets 

sh ij = -igc a (T a ) ik h kj , s 2 h = 0 (69) 

from which the BRST transformation of the Stueckelberg held £“ can be evaluated iteratively, giving 


S £“ = -c a + g -f abc C b p - f amr / mpq ( p £ <7 £ r + Q( g 3j 


(70) 


It is instructive to check here explicitly the BRST invariance of the held C^. For this, it is better to 
employ a matrix notation, namely 


sAf, = —8/j.c + ig[An,c] , sc=-igcc 

sh = —igch , sP = igh^c , 


(71) 


with A n = A < j 1 T a , c = c a T a , £ = £ a T a . From expression (1671) we easily get 

sC^ L = A^h + h\—dfj,c + ig[A^, c])h — igh'A^ ch — h'cd^h + h^d^ch) 

= igh'cA^h — h){d^c)h + igPA^ ch — igPc A^h — igh' A^ch — Pcd^h + P(d fl c)h + Pcd^h 
= 0 , (72) 

which establishes the invariant character of 

The non-local action Sj^ AG , en. (fT3|) . can be now replaced by the equivalent local expression 

J d 4 x - u+ 91 2 e^ ( 4 >f - $?) Cl + r a d,Cl + r 3 ^C 3 ) 


[ S MAg\ loc 


The helds (t",t 3 ) are Lagrange multipliers imposing that the helds (C“,C 3 ) are transverse, i.e. 


d„Cl = 0 


dA = 0 


(73) 


(74) 


As discussed at length in [Tj, see in particular eqs.(A13-A16) of Appendix A, the trasnversality conditions 
(1741) can be solved iteratively for £ as a function of the gauge held A^. In fact, imposing <9 M C“ = <9 M C 3 = 0, 
yields 

d 2 £ = d il A + ig[d il A ll ^]+ig[A ll ,d^\+ g 2 d^A^ + g 2 ^d^A^ + g 2 iA^d^ 

- y 9 “ a A 2 - ~ 9 pM A “ “ 9 A iA “ ~ Y ed “ A “ 


+ *§M 2 £] + 0 (£ 3 ) . 


(75) 
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Solving iteratively, we arrive at 


£ Q2 d/iAn + i q2 


dA, 


dA 

d 2 


+ * 


. 9 


d 2 


Afn dfj, 


dA 

IP 


+ 


* 9_ 
2d 2 


dA 

d 2 


, dA 


+ 0(A 3 


and thus 

where 




(76) 

(77) 


= An - - T^dndA - ig 


dn 


> — “ M Q2 A 1 ”-y Q2 


' , 9A' 

^7/ q2 


— I 


2d 2 


dA, -TpdA 


+ ig 


^A 1 ’ 


. . 9 
+ l 2 


^BA.^BA +0(A 3 ), 


(78) 


so that the starting non-local expression (HD is recovered. 


Let us proceed now with the localization of the BRST operator s 2 , ea. (l46l) . Here, we follow the procedure 
already employed in [93] and we introduce the following auxiliary action 


Sn.UX - 


d A x ^fM aS (C)pff - r/M aS (C)^ - g^pfcl) , 
where are a pair of bosonic fields, while are anti-commuting. 


(79) 


It can be easily checked that the combination ( loc + S aux ) is left invariant by the following local, 

nilpotent, BRST transformations S [ oc : 


SlocA* 

-(DfJ + ge^Afc), 

$locC 

= ge^cPc, s loc c= 9 -e a P 

$locC 

= b h ’ a , s loc c = b h , 

siocb h,a 

= siocb — 0 , 

si oc h lj 

-igc a (T a ) ik h k i , s loc C\ 

Sloe 

= UJ, Sl oc LO“P = 0 , 

Sloc0f 

= +i 2 ^ > si oc 4>^ 

SlocP 

= , siodf = 0 , 


SlocCjj, = 0 


= 0 , 


Su 


= 0 , 


and 
In fact, 


SlocSloc — 0 j Sloe ( [*-’m j 4g] loc ^ au %) — 0 


(80) 

(81) 


sioc([Sl AG ] loc + S aux ) = J d 4 * (<^A4 Q5 (C)a;f - (^ 

+ I d 4 x ( 7 2 PfM a5 (C)ujff - g^ 2 £ af3 cofC 3 ^j 


= 0 


(82) 


As explained in [93], the introduction of the auxiliary action (17911 enables us to localize the BRST operator 
Sry 2 . This follows by looking at the equation of motion of the auxiliary field j3, namely 


SS n 


5/3 


= M(C)/3 — gC 3 = 0 


P = gC 3 


l 


M(C) 


(83) 
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Therefore, from equation (|77|) . one gets 


_ _ 3 1 

siocU = (j) + 7 2 P => s 7 2u) = (f> + g') 2 ( A h ) , 

SlOC ^ Sry2 , 

so that the non-local expression of the operator s 7 , eq. ([^|) , is recovered. 

In summary, the localized actiorjfl 


[^MAg] loc ~ Sym + [ S MAg\ loc + 

+ f d A x (l> h -' a DfAf, - (?M a0 C P + g£ a ^(D^A^)c + b h d ^+ ge a ^A a ^^j , 

(85) 

is left invariant by the local, nilpotent, BRST transformations (18011 . he. 

^ioc [‘S'mAg] loc 0 j SlocSloc 0 . (86) 


Both action (18511 and BRST operator ( 1801 ) reduce to their non-local expressions when the auxiliary 
localizing fields (<^, /3, /5, "0, V’) are integrated out. Moreover, equations ( 1861 ) can be immediately translated 
into Slavnov-Taylor identities useful for studying the algebraic renormalization as well as the cohomology 
of the local BRST operator si oc . We shall post-pone this study to a further detailed analysis. Let us 
end this section by calling attention to the fact that the non-polynomial character of the local action 
[sgioLoc in the Stueckelberg field £ does not jeopardise the use of the powerful tools of the algebraic 
renormalization m , which have been already successfully applied to other non-polynomial actions like, 
for instance, N = 1 Super-Yang-Mills in super-space, chiral Wess-Zumino models in two dimensions as 
well as non-linear sigma models. 


7 Conclusion 


In this work the generalization of the non-perturbative exact nilpotent symmetry of the Gribov- 
Zwanziger action constructed in [I] in the Landau gauge has been generalized to the case of Yang-Mills 
theories quantized in the maximal Abelian gauge, as summarized by eas. (14611 . (1471) . 

It is worth to point out the deep similarity existing between the non-perturbative nilpotent BRST op¬ 
erators of the Landau and maximal Abelian gauges, as one can infer by comparing expressions m and 
m • One notices in fact that the Faddeev-Popov operators of both gauges appear in exactly the same 
way, i.e in the transformation of the auxiliary field oj. The presence of the inverse of the Faddeev-Popov 
operators shows that the non-perturbative BRST operator s 7 2 feels the presence of the Gribov horizon. 
As such, the operator s 7 2 is deeply intertwined with the geometry of the corresponding Gribov regions 
in both gauges. 

The exact non-perturbative BRST symmetry has also been extended to the refined version of the theory, 
ea. (14911 . which takes into account the existence of the dimension two condensates. 


6 We also point out that the Lagrange multipliers (b h , b h,a ) appearing in expression (18511 can be considered as elementary 
fields. This follows by noticing that the field redefinitions in eq. m correspond in fact to field transformations with unit 
Jacobian. 
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The resulting diagonal gluon propagator has been evaluated in D = 4,3, 2 dimensions. While in D = 4, 3 
a decoupling type behaviour has been found, see eas. (|63p . (RJ5L in D = 2 a scaling type behaviour emerges, 
as given in eq. HMD- Once more, we underline the strict analogy existing with the gluon propagator of 
the Landau gauge, which exhibits a similar behavior. This feature suggests a kind of general behavior of 
the gluon propagator in different gauges, as recently advocated in [SH], where a study of the equal-time 
spatial gluon propagator has been performed in the Coulomb gauge, obtaining similar results. 

Finally, we hope that this work will stimulate our colleagues from the lattice community to pursue 
the numerical studies of the diagonal gluon propagator in both D = 4,3, 2 dimensions. 
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